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1. Introduction 

In 1925, Jarmk [3] denned a sequence of convex polygons for use in constructing curves 
containing many lattice points relative to their curvature. Given a positive integer Q, let Vq 
denote the set of all primitive integral vectors in the square of side length 2Q centered at 
the origin, that is, 

Vq = {(q, a) e Z 2 : gcd(g, a) = 1, max{|a|, \q\} < Q}. 

Then the Jarmk polygon Pq is the unique (up to translation) convex polygon whose sides 
are precisely the vectors in Vq. In other words, Pq is the polygon whose vertices can be 
obtained by starting from an arbitrary point in R 2 and adding the vectors in Vq one by one, 
traversing those vectors in a counterclockwise direction. For example, the forty-eight vectors 
in V4, listed in counterclockwise order, are 

(1,0), (4,1), (3,1), (2,1), (3,2), (4,3), (1,1), 

(3,4), (2,3), (1,2), (1,3), (1,4), (0,1), (-1,4), . . . , (1) 

and hence P4 is the tetracontakaioctagon that can be translated to have vertices at 

(-1,0), (0,0), (4,1), (7,2), (9,3), (12,5), (16,8), (17,9), 

(20,13), (22,16), (23,18), (24,21), (25,25), (25,26), (24,30), .... (2) 

These polygons were featured on a recent cover of the Notices of the American Mathematical 
Society in connection with an article of Iosevich [2] and are discussed in further detail in [1, 
Chapter 2] . Figure 1 shows the four sets of vectors V\ through V4 and the four polygons Pi 
through P 4 which they generate. 1 

The polygons Pq have the same eight-fold dihedral symmetry as the unit square E = 
[— 1, l] 2 and, if properly scaled and translated, can be made to pass through the points 
(±1,0) and (0,±1). We denote by Pq these scaled and translated copies of Pq. In a coda 
to [2], Casselman suggests, based on empirical evidence, "that the scaled polygons [{Pq}] 
converge to a somewhat ragged limit curve". The first several scaled Jarmk polygons have 
been superimposed in Figure 2, with a magnified portion shown in the box to the right; the 
darker polygons correspond to larger values of Q. The purpose of this paper is to calculate 
explicitly the limiting curve of the Jarmk polygons. Indeed, we present several variations on 
Jarmk's polygons and calculate the corresponding limit curves, in many cases explicitly and 
in other cases parametrically. 



1991 Mathematics Subject Classification. 52C05 (11H06). 

1 Figure 1 and the boxed portion of Figure 2 are a modification of the cover image for the June/July 2001 
issue of the Notices of the AMS; they were drawn by Bill Casselman in Postscript. 
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Figure 1. Jarnik polygons, right, and their generating sets of vectors, left 




Figure 2. Scaled Jarnik polygons superimposed 

If A is a curve in M 2 and e > 0, let A(e) denote the ^-neighborhood of A, that is, the set of 
all points whose distance to A is less than e. Given a sequence of curves Ai, A^, ■ ■ ■ , we say 
that the curves {Aj} converge to A if for every e > 0, there is some integer j(e) such that 
Aj is contained in A{e) for every j > j(s). Our main result, which we prove in Section 2, is 
the following theorem. 

Theorem 1. Let C be the curve that contains the graph of the equation 

y = \x 2 -l, -§ < x < | 

and that is invariant under rotation by ir/2 around the origin. Then the scaled Jarnik 
polygons {Pq} converge to C . 

The curve C is infinitely different iable everywhere except at the four points (±|,±|), 
where it is only twice differentiable. Although this limiting curve is surprisingly tame, it is 
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Figure 3. The curves C, left, and C l , right 



the case that the local "curvatures" of the scaled Jarmk polygons oscillate rather than tending 
to the corresponding local curvatures of C. We describe this phenomenon in Theorem 7, the 
statement and proof of which appears in Section 6. 

It is interesting to note the relationship between C and the curve Ci, defined as the graph 
of the equation yjl — \x\ + yjl — \y\ = 1 (these two curves are displayed in Figure 3). Indeed, 
if C is rotated by | and then expanded by the factor so that the image again passes 
through the four points (±1, 0) and (0, ±1), then the resulting curve is none other than C\. 
Vershik [5] showed that C\ is the "generic" shape of a convex polygon with lattice point 
vertices, in the following sense: let P be chosen at random uniformly from among all convex 
polygons whose vertices lie in the set (^Z) 2 n E. Then, with probability approaching 1 as 
Q tends to infinity, the polygon P lies within any prescribed open neighborhood of C\. 

That the curves C and C\ differ only up to rotation and scaling suggests that rotating the 
domain demarcating the vectors in Vq might yield interesting results. Very generally, given 
a set S C IR 2 , we may define the set of vectors 

V Q (S) = {(a,b) G Z 2 : gcd(a,6) = 1, ( j, §) G S}. 

We also define the corresponding convex polygons Pq(S) whose sides are precisely the vectors 
in Vq(S), as well as the scaled and translated versions Pq{S) that pass through the four 
points (±1,0) and (0, ±1). If we take S = E, we recover the vectors Vq and polygons Pq 
in Jarmk's original definition above. While these definitions make sense for any set S, it 
seems reasonable in practice to restrict ourselves to sets S that are star-shaped with respect 
to the origin and that are the closures of their interiors; indeed, in this paper we will only 
consider such sets S centered at the origin that in addition have the same eight-fold dihedral 
symmetry as the unit square. We call the Pq(S) generalized Jarmk polygons. 

Let D denote the "unit diamond", namely the square with vertices (±1,0) and (0, ±1). 
The relationship between C and C\ suggests the following theorem, which we establish in 
Section 3. 
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Figure 4. The families O s , left, and B p , right 

Theorem 2. The scaled generalized Jarnik polygons {Pq(D)} converge to C\. 

Put another way, the generalized Jarnik polygons Pq(D) are "generic" in shape, in the sense 
of Vershik's theorem. 

In this paper we also compute the limit curves corresponding to the Pq{S) for two families 
of sets S, both of which were chosen because they interpolate between the unit square E 
and the unit diamond D. 

• For any positive real number S, let Os be the octagon with vertices at (±1, 0) and (0, ±1) 
and the four points (i^, ±y^). These octagons have eight-fold dihedral symmetry, 
and the slopes of the two edges meeting at the vertex (1, 0) are ±5. When 5 = 1 and 
as 5 — > oo, the octagons Os degenerate to the squares D and E, respectively. Figure 
4 shows several of these octagons, with Oin the innermost and Ooo the outermost. In 
Section 4 we calculate the limiting curves of the polygons {Pq(Os)} explicitly, and for 
all values of 5 these limiting curves are comprised of pieces of parabolas. 

• For any positive real number p, let B p be the set {\x\ p + \y\ p < 1}. When p > 1, the 
set B p is simply the closed unit ball in R 2 under the £ p metric. Again, when p = 1 and 
as p — > oo we recover the squares D and E. Figure 4 shows several of these sets, with 
B1/2 the innermost and B^ the outermost. (We remark that the boundary of Bi/ 2 is 
also closely related to Vershik's curve C\.) In some cases, we can explicitly compute the 
limiting curves of the Pq(B p ), and in these cases the limiting curves are again piecewise 
algebraic. In all cases, we obtain a parametric representation of the limiting curves and 
suspect that they are not in general piecewise algebraic. This family of examples is 
investigated in Section 5. 

We note in passing that one can consider the problem of generalizing Vershik's theorem to 
domains other than the unit square E. Given a set S G R 2 , choose a polygon P at random 
uniformly from among all convex polygons whose vertices lie in the set (^Z) 2 fl S. Is there a 
curve V(S) such that, with probability approaching 1 as Q tends to infinity, the polygon P 
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lies within any prescribed open neighborhood of V(S)? Vershik's result is that V(E) = C\, 
it seems likely that V(D) is the curve C scaled by a factor of |, so that it passes through the 
points (±|, ±|). In general, it would be interesting to search for a connection between these 
generalized "Vershik curves" V(S) and the limiting curves of familes of generalized Jarnik 
polygons {Pq(S')} for appropriate sets S and S'. 



2. The Original Jarnik Polygons 

Because of the eight-fold symmetry of the Jarnik polygons, we need only consider the 
portion of Pq starting from the edge corresponding to the vector (1,0) and ending with 
the edge corresponding to the vector (1, 1); we call this eighth-portion the fundamental arc 
of Pq, as the entire polygon Pq is generated from the fundamental arc under the action 
of the dihedral group of order eight. For example, the fundamental arc of P 4 consists of 
the seven edges defined by the first eight vertices in equation (2). Given any curve C with 
eight-fold dihedral symmetry about the origin (including the scaled Jarnik polygons Pq and 
their generalizations), we shall also refer to the eighth-portion of the curve lying in the wedge 
{(x, y) : x > 0, y < —x} as the fundamental arc of C. 

For any real number A G [0, 1] we define Vq(X) to be the set of all vectors in Vq with positive 
coordinates and slope not exceeding A. The sum of all the vectors in Vq(A) corresponds to 
a particular vertex on the fundamental arc of Pq. If we translate Pq so that the right-hand 
endpoint of the edge corresponding to the vector (1,0) is at the origin, as in (2), then the 
coordinates (X(Q, A), Y(Q, A)) of this vertex are given by the formulas 

X(Q,A) = E<? E 1 and Y(Q,\) = J2 E a - 

q<Q a<Xq q<Q a<\q 

gcd(<?,a)=l gcd(cj,a)=l 

For example, we see from equation (1) that V A (^) consists of the three vectors (4,1), (3,1), 
and (2,1), and hence (X(4, ^),Y(4, ^=)) is the vertex (9,3) of P 4 . 

The following asymptotic evaluation of (X(Q, A) and Y(Q, A)) is the key to our calculation. 



Lemma 3. We have 



X(Q, A) = + 0(Q 2 log Q) and Y(Q, A) = + 0(Q 2 XogQ) 



uniformly for Q > 2 and < A < 1 . 

Proof: Recall the definition of the Mobius mu-function 
1, if n = 1, 

/i(n) = ^ (— l) r , if n — P1P2 . . .p r where the pt are distinct primes, 
0, if the square of any prime divides n. 

The well-known Mobius inversion formula is based on the characteristic property of n 

YvlW = { 1 ' ifn = 1 ' (3) 

4^ |0, if n> 1. K J 
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It is also well known that Y2q>i M?)/? 2 = VC(2) = 6/7r 2 , where ( denotes the Riemann 
zeta-function. We shall use the truncated version of this identity 

which follows easily by a trivial estimation of the tail by J2 q >q V<? 2 - 
Property (3) allows us to write 

X(Q,A) = £gW £ M<0) E 1 - 

d|g d\a 

Writing a = bd and g = rd, we have 

X(Q,A) = £M^) £ dr£l 

d<<2 r<Q/(2 b<Ar 

= X>(d) £ r(Ar + 0(l)) 

rf<Q r<Q/d 

Jl , /O 1 \ ^ j \ ^ _ \ 

(5) 



= A £>(<*) £ r 2 + £rf £ r 

(2<Q r<Q/d ^d<Q r<Q/d 



AE« d )(i(§) 3 + o((f) 2 )) + o(E4) 2 ) 



Equation (4) now allows us to conclude 

X(Q,X) = ^(- 2 +0(h)+0(Q 2 logQ) = ^ + 0(Q 2 log Q) 

as claimed. 

In the same way we see that 



mA) = ££«( £ H(d)) 



q<Q a<Xq d\gcd(q,a) 

= £/^)££° 

d<Q q<Q a<Xq 

d\q d\a 

= £Md) £ J2 db 

d<Q r<Q/db<Xr 

= £>(<*) £ (^(Ar) 2 + 0(Ar)) 

d<<2 r<Q/d ^d<Q r<Q/d 
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At this point, a direct comparison to the middle line of equation (5) yields 

Y(Q,\) = ^f + 0(Q 2 \ogQ) 

as claimed. □ 

We can now prove Theorem 1. Define R(Q) = X(Q, 1) + Y(Q, 1) - 1/2, so that R(Q) = 
3Q 3 /ir 2 + 0(Q 2 \ogQ) by Lemma 3. If we translate Pq so that the midpoint of the edge 
corresponding to the vector (1, 0) is at the point (0, —R(Q)), then the center of Pq will be at 
the origin due to the symmetries of Pq; we then obtain Pq by scaling by the factor 1/R(Q). 
If (X(Q, A), Y(Q, A)) is the vertex of Pq corresponding to the vertex (X(Q, A), Y(Q, A)) of 
P Q , then 

using Lemma 3. In particular, when Q is large enough, the vertices on the fundamental arc 
of Pq lie within e/2 (say) of the arc parametrized by (2 A/3, A 2 /3 — 1) with < A < 1. This 
parametric curve is precisely the arc of the parabola y = 3x 2 /A — 1 from x = to x = 2/3, 
which is an eighth-portion of the curve C. Moreover, the lengths of the edges of Pq are 
0(1/Q 2 ), and so every point on the fundamental arc of Pq lies within e of C when Q is large 
enough. Finally, because of the symmetries of C and the Pq, we see that the entire polygon 
Pq lies within an e-neighborhood of C when Q is large enough. This establishes Theorem 1. 



3. Polygons Defined by the Unit Diamond 

We begin by looking at the derivation of Lemma 3 from another viewpoint. For any real 
number < A < 1, define E(X) = {(x,y) £E:x>0,0<y< Xx}. The inner double sum 
in the first line of equation (5) is written as a sum over lattice points in a large wedge, but 
we may reinterpret it as a sum over (^Z,) 2 fl -E'(A) by writing 

E"M£*£i = <? 3 E^((h) 2 E § E 

d<Q r<Q/d b<Xr d<Q ^ V r : 0<dr/Q<l ^ b : 0<db/Q<Xdr/Q ' 

Notice that the quantity in parentheses is a Riemann sum approximating the integral 




x dx dy, 
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and in fact (since the integrand x has bounded first derivatives) the error in making this 
approximation will be proportional to the mesh size, which is 0(d/Q). Therefore 



x(Q,\) = Q 3 J2^r( [[ xdxdy + od) 



" C(2) J J EW 

This is in agreement with Lemma 3, as ((2) = tt 2 /6 and 



II xdxdy + 0{Q 2 log Q). 

J J E(X) 



xdxdy= I (I xdxjdy = —. 

E(X) Jo \Jy/X J 3 

Similar remarks apply to Y(Q, A). 

We shall use a similar approach for the generalized Jarmk polygons Pq(S). For any real 
number A G [0,1] we define Vq(S, A) to be the set of all vectors in Vq(S) with positive 
coordinates and slope not exceeding A. The sum of all the vectors in Vq(S, A) corresponds 
to a particular vertex on the fundamental arc of Pq(S). If we translate Pq(S) so that the 
right-hand endpoint of the edge corresponding to the vector (1, 0) is at the origin, then the 
coordinates (Xs(Q, A), Y$(Q, A)) of this vertex are given by the formulas 

X 5 (Q,A) = ]Tg Yl 1 and Y(Q,\) = J2 E a - 

q<Q (q,a)€V Q (S) q<Q (q,a)eV Q (S) 

a<\q a<Xq 

If we define 5(A) = {(x, y) <E S : x > 0, < y < \x}, then the same argument as above 
allows us to conclude that 

X S (QA) ~ J J xdxdy and Y S (Q, A) ~ Jj ^ y dx dy. (6) 

This depends of course on S being a "reasonable" set. For the sets S we shall consider, the 
asymptotic formulas (6) do in fact hold, with error terms that are 0(Q 2 \ogQ). 
We also use the definitions R S {Q) = X S {Q, 1) + Y S (Q, 1) - 1/2 and 

~ X s (Q,\) + l/2 ~ Y S (Q,\)-R S (Q) 
X5(g ' A)= Rs(Q) ' Ys ^ X) = lUQ) ' 

so that (X S (Q,\),Y S (Q,\)) will be the coordinates of the corresponding vertex on the 
fundamental arc of the scaled generalized Jarmk polygon Pq(S), the translation and scaling 
chosen so that the center of Pq(S) is the origin and the points (±1,0) and (0, ±1) are 
midpoints of edges of Pq(S). 

We now implement this approach with S equaling the unit diamond D to prove Theorem 2. 
Again, due to the symmetries of C\ and the Pq(D), it suffices to show that the fundamental 
arcs of the Pq{D) tend to the fundamental arc of C±. Given < A < 1, the set D(X) is the 
same as {(x, y) : x > 0, x + y < 1, y < \x} = {(x, y) : < y < A/(l + A), y/X < x < 1 — y}. 
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Therefore 

JJ Dm ydxdy= l y [L dx r = wrw 

Using these evaluations in equation (6), we see that 

^,A)~ Q ;y + ^ , y D ( Q ,x)- Q3x2 



7T 2 (1 + A) 2 ' ' 7T 2 (1 + A) 2 

as C(2) = vr 2 /6. This implies that R D (Q) ~ Q 3 /ir 2 , and so 

X D (Q,X) + l/2 A(2 + A) 



X D {Q,\) 
Y D (Q,X) 



R D (Q) (1 + A) 2 

Y D (Q,X)-R D (Q) -(2X + 1) 



(7) 



R D (Q) (1 + A) 2 • 

If we set x = jpf^p and y = -j^^r , it is easy to check that a/1 — fx] + a/1 — \y\ = 1, and 
hence the curve parametrized by ( (i + ~^2 , '(f+xf ) with < A < 1 is precisely the fundamental 
arc of C\. This establishes Theorem 2. 

4. Polygons Defined by Octagons 

Recall that for any positive real number 5, we defined 0$ to be the octagon with vertices 
at (±1,0) and (0, ±1) and the four points (±^, ± j-f^)- We can use the same strategy to 
calculate the limiting curve of the generalized Jarmk polygons generated from the sets 0$. 
Define C$ to be the curve with eight-fold dihedral symmetry whose fundamental arc is 

{(*, V) ■■ < x < |f±i, 45(1 + 5) 2 (y + 1) = (1 + 35)(5x + y + l) 2 }. (8) 
This arc is part of a parabola whose axis of symmetry has slope —5 and whose vertex 
is ( (1+3^(1^2)2 ; ^ 1 |i+3/)^i+^2)2 3 ' 5 ' )) as it turns out. The endpoints of this parabolic arc 
are (0,-1) and (||ty, — ffrj)- 111 particular, each C$ is piecewise algebraic, and one can 
check that each Cs is twice differentiable at the eight symmetry points (±1, 0), (0, ±1), and 
(±||^Y, ±||±y). Figure 5 shows several of these curves, with C1/4 the outermost and the 
innerermost; the points (±|,±|) and (±|,±|), which lie on and Ci, respectively, are 
also indicated. 

Theorem 4. For every positive real number 5 , the scaled generalized J arnik polygons {Pq{0 5)} 
converge to Cs- 

Note that when 5 = 1, the octagon 0\ degenerates to the unit diamond D; in this case, 
the equation of the parabola in equation (8) is equivalent to4(y + l) = (x + y + l) 2 , which 
turns out to be another way to define the fundamental arc of the curve C\. Therefore the 
notation Cs is consistent with our earlier definition of Ci, and Theorem 4 is consistent with 
Theorem 2. Also, as 5 tends to infinity, the octagons 0$ converge to the unit square E, while 
the equation of the parabola in equation (8) tends to 4(y + 1) = 3a; 2 , which is the equation 
defining the fundamental arc of C . Therefore Theorem 4 has Theorem 1 as a limiting case 
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as well. It can also be checked that C$ tends towards the boundary of the unit square E as 
5 decreases to zero. 

We prove Theorem 4 using the same approach as the proof of Theorem 2 in the previous 
section, showing the important steps while supressing the algebraic details of the computa- 
tions. Given < A < 1, we have by equation (6) 

x o 5 (QA)~-§7^[[ xdxdy 

U 2 J J J Os(X) 
q3 r 8X/(8+X) / r l-y/S 



C(2) 



,6W+\) / rl-y/6 X ^ Q 3 5\(25 + A) 

Jo \L xdx ) dy -^ (5 + xy 



Yo„(Q,\) ~ / / ydxdy 



OsW 

Q 3 f S ^+X) / rl-y/S X Q3 5 2 X 2 

/ y[ dx\dy 

JO \ Jy/X J 



This implies that Rq s (Q) ~ ^r^j^r, and so 



C(2)7 *\J y/ x J * ^(s + xy 

_ . ) o n r\ ar\ 

{5+iy 



~ X 0s (Q, A) + 1/2 A(2<5 + A)Q5+1) 2 

Xoa(g ' A) - i2o a (g) FWTI) 
~ r ,(g,A)-^(g) <5a 2 (^ + i) 2 

If we set x = X (s+\y(38+i) an< ^ V = (<5+a)^(3^+i) ~~ ^ one can cnec k that { x iV) satisfies the 
polynomial relation (8), and hence the curve parametrized by ( X (s+xp(3S+i) ' (<5+aP(3^+i) ~~ -0 
with < A < 1 is precisely the fundamental arc of C$. This establishes Theorem 4. 
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5. Polygons Defined by Unit Balls 

Recall that for any positive real number p, we defined B p to be the set {\x\ p + \y\ p < 1}, 
which we refer to as the "unit £ p -ball" (an abuse of notation when p < 1). We also need 
the standard notation B(a,b) for the Euler beta function B(a,b) = ^t a ~ l {\ — i) b ~ l dt = 
T(a)r(6) /r(a + 6) as well as its relatives, the incomplete beta function B z (a, b) = J* t a_1 (l — 
t)^ 1 dt and the regularized incomplete beta function I z (a, b) = B z (a, b)/B(a, b). 

Let p be a positive real number which we regard as fixed. When < A < 1, set /i — /i(A) = 
jr^p, and define C' p to be the curve with eight-fold dihedral symmetry whose fundamental 
arc is given parametrically by 

1 + 1)- \ B n L i^i + l)- B( i g 1), 0<A<1. (9) 

Figure 5 shows several of these curves, with C(, 3 the outermost and C'^ the innerermost; 
the points (±|, ±|) and (±|, ±|), which lie on and C[, respectively, are also indicated. 

Theorem 5. For every positive real number p, the scaled generalized Jarnik polygons {Pq(B p )} 
converge to C' p . 

Although the curves C' p are in general rather inscrutable, it can be shown that they are all 
twice differentiable at the eight points of symmetry and infinitely differentiable everywhere 
else. In the special case p — 1, the ball B\ is exactly the unit diamond D. The parametric 
representation (9) of the fundamental arc of C[ reduces to ( (i+tp ' "(i+a)^ ) > wn i cn by equation 
(7) is the parametric representation of the fundamental arc of C\. (Again in this section, 
we supress the details of many of our calculations.) Therefore Theorem 5 is consistent with 
Theorem 2. 

It can also be shown that as p tends to infinity, the parametric representation (9) of 
C' p approaches (2 A/3, A 2 /3 — 1), which is the parametrization of the fundamental arc of C. 
Therefore Theorem 5 has Theorem 1 limiting; well. Again, it can be checked 

that C' p tends towards the boundary of the unit square E as p decreases to zero. 

The case p = 2 is also special, since the domain B 2 (the unit disk) has complete rotational 
symmetry. Indeed, the parametric representation (9) of the fundamental arc of C' 2 reduces to 
(yT+A 5 ' ~ v / ITa t ) ' wn i cn i s the fundamental arc of the unit circle. Thus the limiting curve of 
the generalized Jarnik polygons formed from the vectors in Vq(B2) is simply the unit circle, 
not surprisingly. 

When p is the reciprocal of a positive integer, the regularized incomplete beta functions 
are simply indefinite integrals of polynomials. Therefore in these cases, the parametric 
representation (9) of the fundamental arc of C' p can be written as rational functions of X p . 
In particular, these particular curves C are piecewise algebraic, and in principal one can 
calculate the algebraic equation defining the fundamental arc. For example, when p = 1/2 
the parametric representation (9) of the fundamental arc of C[^ 2 reduces to 

/ Aqo + IOA^ + SA + A 3 / 2 ) -(l + SA^ + lOA + lOA 3 / 2 ) ^ 

V (1 + AV2)5 ' (l + A 1 /2)5 J' °< A < 1 ' 
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and the coordinates (x, y) of this parametrization satisfy the irreducible polynomial relation 

-45253 + 86140a; - 37030a; 2 - 3220a; 3 - 765a; 4 + 128a; 5 - 86140?/ + 169060x2/ 

- 80340a; 2 y - lM0x 3 y - 6A0x 4 y - 37030y 2 + 80340a;?/ 2 - 44590a; 2 ?/ 2 + 1280a; 3 2/ 2 

+ 3220y 3 - IMOxy 3 - 1280a;V - 765y 4 + 6A0xy 4 - 128y 5 = 0. 

We establish Theorem 5 using our now familiar technique. For any < A < 1, equation 
(6) gives 




This implies that 

r b p (Q) ~ §T)ik B (h §) = mh B ^ 1 + §) = §kp B (h 1 + p)' 

and so after much calculation we see that 

x Bp (Q, A) ~ i + 2(1 - AP p PjB( i ? I) 

Y Bp (Q, A) ~ -[h-,( l p , l + j)- 2 (i + AP)3/pB(i,|))' 

which is exactly the parametric definition (9) of the fundamental arc of C' p . This establishes 
Theorem 5. 

6. Local Curvatures 

Given a vertex v of any polygon P, we define the radius of curvature of P at v to be the 
radius of the circle passing through v and its two neighbors. We quantify the local curvatures 
of the Jarmk polygons (as originally defined) in the following way. To each irrational number 
< A < 1, we associate the unique vertex vq(X) on the fundamental arc of Pq such that 
A lies between the slopes of the two edges adjacent to vq(X). We then define tq(X) to be 
the radius of curvature of Pq at v(Q,X). This description is not well-defined for rational 
numbers A, but we can speak of tq(X + ) and tq(X~). For example, from equation (2) we see 
that r Q(-^) is the radius of the circle passing through the points (7,2), (9,3), and (12,5), 

which turns out to be ^1105/2. We also have r Q (\ + ) = ^1105/2 but r Q (\~) = 5^29/2. 

After scaling the Jarmk polygons, the radius of curvature tq(X) at the corresponding 
vertex of Pq is simply tq{X)/ R{Q). It would be tidy if, as Q grew large, the local radii of 
curvature tq(X) would converge to the radius of curvature of the limiting curve C at the 
corresponding point (2A/3,A 2 /3 — 1), which turns out to be |(1 + A 2 ) 3 / 2 . However, not 
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10 10 2 10 3 10 4 10 5 10 6 10 10 2 10 3 10 4 10 s 10 6 1 7 1 8 

Figure 6. The local radii of curvature r Q (A) for A = 4=, left, and A = e — 2, right 



only does {^q(A)} never converge to |(1 + A 2 ) 3 / 2 , but in fact {^q(A)} fails to converge at 
all for most A, and the manner in which it fails to converge depends upon the diophantine 
approximation properties of A. In Figure 6 we have plotted these local radii of curvature 
r<g(A) as functions of Q (represented on the horizontal axis in logarithmic scale) for two 
interesting examples of irrational numbers, A = -j= and A = e — 2, with the horizontal 

dashed line indicating the value |(1 + A 2 ) 3 / 2 in each case. 

We recall some notation and standard facts from the theory of diophantine approximation 
and ontinued fractions. The Farey fractions of order Q are defined to be the rational numbers 
in [0, 1] with denominator not exceeding Q, listed in increasing order. For example, the Farey 
fractions of order 4 are {f>|>|>§>§>f>Y}- There is a one-to-one correspondence between the 
Farey fractions - of order Q and the vectors (q, a) in the "fundamental arc" of Vq, as we see 
from equation (1) when Q = 4. If ^ and ^ are consecutive Farey fractions, it is known (see 
[4, Section 6.1]) that a^q\ — a\q2 = 1. 

Let A be an irrational number in (0,1) with continued fraction expansion [0; bi, 62, • • • ] 
where the b n are positive integers. Define sequences {h n }, {k n } of positive integers by 

ho =1, h x = 0, h n+ i = b n h n + h n ^ (n > 1) 
k = 0, h = 1, k n+ i = b n k n + £; n _i (n > 1). 

The {h n /k n } are the convergents to A. It is immediate that for any real number Q > k,2 — bi, 
there is a unique index n > 2 and a unique integer 1 < j ' < b n such that jk n + k n _i < Q < 
(j + l)k n + k n -i. It is known that in the set of Farey fractions of order Q, the number A lies 
between two fractions whose denominators are k n and jk n + fc n _i in this notation. (See [4, 
Section 7.5, Problem 5]. Fractions of the form (jh n + h n -i)/(jk n + fc n _i) with 1 < j < b n 
are called the secondary convergents to A.) 

For example, if A = A= — [0; 1, 1, 2, 1, 2, 1, 2, . . .], then the sequence of convergents is 

{ i> T' \i I' I' fl> M ' • • • }■ Setting Q = 15, for instance, we have A; 4 = 5, k 5 = 7, j = 1 < 
2 = 65, and 1-7 + 5 < Q < 2-7 + 5. Hence the largest (respectively, smallest) rational 
number with denominator bounded by 15 that is less than (respectively, greater than) 

is I (respectively, = ^), and therefore the two edges in the fundamental arc of P15 

between whose slopes 4= lies correspond to the consecutive vectors (7,4) and (12,7) of Vi 5 . 



constant 5 > such that the inequality |A — || < has no solutions. The set of badly 
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An irrational number < A < 1 is badly approximable if the partial quotients bj in the 
continued fraction expansion A = [0; bi, b 2 , • • • ] are bounded, or equivalently if there is a 

q I < 

approximable irrationals has Lebesgue measure zero. 
Lemma 6. For any irrational number < A < 1, we have liminf(/c n /A; n+1 ) < (v^5 — l)/2. 

Proof. Whenever b n > 2 we have k n /k n+ \ = k n /(b n k n + k n -i) < l/b n < 1/2. Therefore 
if infinitely many of the b n > 2, then liminf (k n /k n+ 1) < 1/2. Otherwise, b n — 1 for 
n sufficiently large, so the /c n eventually satisfy k n+ i = k n + k n -\. All solutions to this 
recurrence in positive numbers satisfy k n ~ c((v / 5 + l)/2) n for some constant c, hence 
lim(W*»+i) = ((v^ + l)^)- 1 = (v/5 - l)/2. □ 
We can now describe the limiting behavior of the local radii of curvature tq{X). 

Theorem 7. Le£ < A < 1. 

a. If X is rational, then \im.Q^ OCl fQ(X + ) = limQ^oo fg(A _ ) = 0. 

b. If X is irrational, then limsupg^^f^A) lies in the interval 

[^(1 + A 2 ) 3/2 ,^(1 + A 2 ) 3/2 ]- 

c. liminf q^qo tq(X) > if and only if X is a badly approximable irrational number. 
In particular, for almost all X, we have limsupQ^^ fg(A) > but liminfg^oo fg(A) = 0. 

Proof. A straightforward calculation shows that the radius r of the circle passing through 
the three points (x — x±, y — y±), (x, y), and (x + X2, y + y 2 ) satisfies 

r 2 = \{yl + x\)(yl + 4) {{yi + yif + (a* + x 2 ) 2 ) (y 2 x 1 - yi x 2 )- 2 . (10) 

To calculate tq(X), we take (xi, y±) = (qi, ai) and (x 2 , y 2 ) = (q 2 , a 2 ), where A lies between ^ 
and ^ in the Farey fractions of order Q, and (x,y) = {X(Q, A), Y(Q, A)). Since ^ and ^ 
are consecutive Farey fractions, we know that a 2 q\ — a±q 2 = 1, and hence the formula (10) 
simplifies to 

r Q {Xf = \{a\ + ql)(a 2 2 + q\) ((a, + a 2 f + ( qi + q 2 f). (11) 

Now ^ and ^ are both approximately A, so substituting a\ ~ Agi and a 2 ~ Ag 2 into equation 
(11) and simplifying yields 

rQ{X)^\q l q 2 {q l + q 2 ){l + X 2 fl 2 . 

(We record only the main terms for the sake of simplicity. Even a crude estimate such as 
\j- — A I = o(^) would suffice for our purposes.) Therefore 

- m r Q (X) gig 2 (gi + g 2 )7r 2 (l + A 2 ) 3 / 2 

r ° (A) - mg) 6 — • (12) 

If A = - q is a rational number, then in using the expression (12) to calculate fg(A + ), 
we would have q\ = q for all Q > q. In particular, since q 2 < Q, the numerator grows 
only quadratically with Q, and hence limQ^oo fg(A + ) = 0. The same argument shows that 
limg^oo fg(A~) = 0, which proves part (a) of the theorem. 

Now suppose that A = [0; bi, b 2 , . . . } is irrational. Then there are unique positive integers 
n and j < b n such that jk n + k n ^i < Q < (j + l)k n + /c n _i, where the k n are the denominators 
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of the convergents to A, in which case qi = k n and g 2 = jk n + k n _\. Thus the expression (12) 
becomes 

~ / % fc„(jfc n + fc„_i)((j + 1)A;„ + /c n _i) tt 2 (1 + A 2 ) 3 / 2 

rg(A) ^ g . (13) 

In calculating the lim sup of this expression, we should take Q as small as possible, that is, 
Q = jk n + k n _i. If we define r n = k n -i/k n , then the expression (13) simplifies to 

U(j + l)k n + fcn-l) 7T 2 (1 + A 2 ) 3 / 2 j + 1 + T n 7T 2 (1 + \ 2 f 2 



(jk n + k n -i) 2 6 (j + r n ) 2 6 



The expression = — I- p is a decreasing function of both j and r n , so in 

calculating the lim sup it is best to take j = 1 and r n as small as possible. Therefore 

2 + liminf r n tt 2 (1 + A 2 ) 3 / 2 

lim sup rQ\X) 



1 + liminf r ra ) 2 6 

(l+(V5-l)/2)2 ' (1+0) 2 



and by Lemma 6 the first fraction lies in the interval [ (i+^^y^ > (l+op ] = [1>2]. This 
establishes part (b) of the theorem. 

Similarly, in calculating the lim sup of the expression (13), we should take Q as large as 
possible, that is, Q = ((j + l)k n + k n _i)~ , in which case (13) simplifies to 

- m UA + U tt 2 (1 + A 2 ) 3 / 2 = j + r n tt 2 (1 + A 2 ) 3 / 2 
r ° U ~ (0' + l)*n + *»-i) 2 6 (j + l + r„) 2 6 

The expression ^ r^" p always lies between J^) 2 an< ^ (j+ip < J' wn i cn are decreasing 
functions of j, so in calculating the lim inf it is best to take j as large as possible, that is, 
j = b n . If A is badly approximable, so that b n < B for all n, then liminf tq(A) > j§^i > 0; 
if A is not badly approximable, then the b n are unbounded above and hence 

lim inf f Q (A) < liminf- y — — '- — = 0. 

bn 6 

This establishes part (c) of the theorem. □ 

The two examples in Figure 6 illustrate the two possibilities in part (c) of the theo- 
rem. As noted before, the continued fraction expansion of is [0; 1,1,2,1,2,1,2,...], and 

so in particular ^= is badly approximable since the partial quotients are bounded above 
by 2. We can see the repeating groups of a single curve followed by a pair of curves in the 
left-hand graph; in particular, the near-periodicity of the graph implies that the values of 
fg(-i=) are bounded below. On the other hand, the continued fraction expansion of e — 2 is 

[0; 1, 2, 1, 1, 4, 1, 1, 6, 1, 1,8,...], and in particular e — 2 is not badly approximable since the 
partial quotients are unbounded. In the right-hand graph we can see the influence of these 
partial quotients (the last full group contains two single curves and a group of 14 curves, 
corresponding to the string 1, 1, 14 in the continued fraction), and in particular that the lim 
inf of the values of fg(e — 2) is zero. 

We remark that the possible values for liminf(/c n /A; ri+1 ) in Lemma 6 are closely related 
to the Markov spectrum. Therefore the possible values for limsupQ^^ fg(A), as well as the 
possible values for liminfg^oo rg(A) for badly approximable irrationals A, are also related to 
the Markov spectrum. 
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